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ABSTRACT 

Super-Eddington accretion discs with 3 and 15 Me around black holes with 
mass IOMq are examined by two-dimensional radiation hydrodynamical cal- 
culations extending from the inner disc edge to 5 x lO^rg and lasting up to 
~ lO^rg/c. The dominant radiation-pressure force in the inner region of the 
disc accelerates the gas vertically to the disc plane, and jets with 0.2 - 0.4c are 
formed along the rotational axis. In the case of the lower accretion rate, the 
initially anisotropic high- velocity jet expands outward and becomes gradually 
isotropic flow in the distant region. The mass-outflow rate from the outer 
boundary is as large as ~ 10^^ - 10^"^ g s~^, but it is variable and intermittent 
with time; that is, the outflow switches occasionally to inflow in the distant 
region. The luminosity also varies as ~ 10^° -10^^ erg s^^ on a long time-scale. 
On the other hand, the jet in the case of the higher accretion rate maintains 
its initial anisotropic shape even after it goes far away. The mass-outflow rate 
and the luminosity attain to steady values of 3 x 10^^ g s"^ and 1.3 x 10*° erg 
s^^, respectively. In accordance with the local analysis of the slim accretion 
disc model, the disc is thermally unstable in the case of 3 A/e but stable in 
the case of 15 A/e- The super-Eddington model with 15 A/e is promising to 
explain a small coUimation degree of the jet and a large mass-outflow rate 
observed in the X-ray source SS 433. 
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1 INTRODUCTION 



Disc accretion is an essential process for such dynamic phenomen a as energetic X- r ay sou rces, 



active galactic nucle i , and protostars. Since the early works by iPringle fc ReesI ()l972n and 
Shakura fc SunvaevI (|l973l ). a great number of papers have been devoted to studies of the 
disc accretion onto gravitating objects, and the standard accretion disc models with geo- 
metrically thin disc have proven to be particularly successful in applications to cataclysmic 
variables. However, for very luminous accretion discs whose luminosity exceeds the Edding- 
ton luminosity, we need another disc model instead of the standard model because the 
disc is expected to have geometrically thick structure and powerful outflows from the disc. 
It is well known that the luminous accretion discs, where the radiation-pressure and the 
electron-scattering are dominant, are thermally unstable if the usual a-model for the vis- 
cosity is used. For the lu minous disc, the slim accretion disk model has been proposed by 



Abramowicz et al, 



()l988f ). where strong advective cooling depresses the thermal instability. 
The slim disc model seems to be reasonable for the accretion discs at the super-Eddington 
luminosity, but we need to reconfirm it from the point of view of a two-dimensional simula- 
tion. The super-Eddington accretion discs were generally expecte d to possess vortex funnels 



and radiation-pressure driven i( 


Lvnden-Bell 


1978; 


Fukue 


11982) 



ets from geometrically thick discs (|Shakura &: Sunvaev 



1973; 



4). A two-dimensional hydrodynamical c alculation of a 



Eddington accretion disc around a black hole was first examined by lEggum et al 



super 



(ma 



19881 ). and was discussed regarding SS 433. Although they showed a relativistic jet forma- 
tion just outside a conical photosphere of the accretion disc and a collimation angle of ~ 30° 
of the jet, the jet's mass fiux is too small and the collimation angle is much larger than the 
one observed for SS 433. 



Blandford fc BegelmanI ()l999[ ) investigated an adiabatic infiow-outfi ow solution of the 
accret ion disc, where the super-Eddington fiow leads to a powerful wind. iKing fc Begelman 
( I999I ) discussed the super-Eddington mass transfer in binary systems and how to accrete in 



such super-Eddington accretion discs. Rec ent theoretical studies and two-dimensional sim- 



ulations of the black holes accretion fiows (llguraenshchev &: Abramowicz 



1999 : 



Naravan et al .120001 



Ouataert fc Gruzinovll200C ; 



Stone fc Pring1ell200ll 



1999; 



Stone et al 



Abramowicz fc Igumenshchev 



20011 ) showed a new development in the advection-dominated accretion fiows (ADAFs) and 
the convection-dominated accretion fiows (CDAFs) which may be deeply related to the 
super-Eddington accretion fiows. Furthermore, from a point of view completely different 
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from the a-viscosity models for accretion discs, recently developed rigorous works of global 
three-dimensional magnetohydrodynamic f MHD) simulations of b lack hole accretion flows 



should be noticed faawley fc Balbui 



2mm 



Ralbus fc Hawlev 



20031 ) ■ They show that, with- 



out assumption of the phenomenological a-yiscosity, magnetic fields play the most important 
role in the angular momentum transport; that is, accretion is driven by turbulent stresses 
generated selfconsistently by the magnetorotational instability. The results reveal a three 
cor nponent structure: a warm Keplerian disc, a highly magnetized corona, and axial jet. 



Okudal ()2002r ) examined the super-Eddington black-hole model for SS 433 based on two- 



dimensional hydrodynamical calculations coupled with radiation transport. Although, in 
this model, a high-velocity jet with 0.2 - 0.4c is formed along the rotational axis, a small 
collimation angle of ~ 0.1 radian of the jets and a sufficient mass-outflow rate comparable 
to ~ 10^^ g s^^ could not be obtained. Additionally, it was also suggested th at this result 



might be a transient behavior during a long disc evolution (jTeresi et al 



20041 ). Keeping in 



mind these facts, we examine further super-Eddington accretion discs and jets around the 
black hole extending from the inner disc edge to 5 x lO^rg. 



2 MODEL EQUATIONS 



A set of relevant equations consists of six partial differential equations for density, mo- 
mentum, and thermal and radiation energy. These equations include the full viscous stress 
tensor, heating and cooling of the gas, and radiation trans port. The radiation transport i s 
treated in the gray, fiux-limited diffusion approximation (jLevermore fc Pomraningill981l ). 
We use spherical polar coordinates {r,(,f), where r is the radial distance, ( is the polar 
angle measured from the equatorial plane of the disc, and if is the azimuthal angle. The gas 
flow is assumed to be axisymmetric with respect to Z-axis {d/d(f = 0) and the equatorial 
plane . In this coordinate system, the basic equations for mass, rnomentum, gas energy, and 
radiation energy are written in the following conservative form |K1evlll989l ): 
dp 
di 



div{pv) = 0, 



d{pv) 

d{prw) 
dt 



+ div{pvv) = p 
div{prwv) = 



(r - rg)2 



dp 1 

- T- + fr + divSr + -Sr 

or r 



dp 



pv^tanC - ^ + div(r5^) + S'.^^tanC + /(, 



d{prcos(Vy,) 
di 



-|- div(prcosCf<^i') = div(rcosC'S'^), 



(1) 
(2) 

(3) 
(4) 
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dps 



dt 
and 



+ <l\v{pev) = —p divv + $ — A, 



+ divFo + div{vEo + v ■ Pq) = A - p 



K + cr 



-V ■ Fo, 



(5) 



(6) 



dt c 
where p is the density, v = {v, w, v^) are the three velocity components, G is the gravitational 

constant, is the central mass, p is the gas pressure, e is the specific internal energy of 

the gas, Eq is the radiation energy density per unit volume, and Pq is the radiative stress 

tensor. The subscript "0" denotes the value in the com oving frame and that the equations 



are correct to the first order of i;/c (Katoetal 



1998^ . We adopt the pseudo-Newtonian 



potential (jPaczvhsky fc Wiitalll98fll ) in equation (2), where rg is the Schwarzschild radius. 



The force density /r = (fr, f() exerted by the radiation field is given by 

K + a 



/r = P- 



-Fo, 



(7) 



where k and a denote the absorption and scattering coefficients and Fq is the radiative fiux 
in the comoving fr ame. For the opac it ies k we use polynomial fits to the Rosseland mean 
Lin &: Papaloizoul (|l985| ). S = {Sr, S(^, S^) denote the viscous stress 



opacities given by 



tensor which includes all components due to the three velocity components (|Okuda et al 



19971 ). $ = (5 V)f is the viscous dissipation rate per unit mass. 



The quantity A describes the cooling and heating of the gas, 

A = pcK{S^- Eo), (8) 

where S"* is the source function and c is the speed of light. For this source function, we assume 
local thermal equilibrium = aT'^, where T is the gas temperature and a is the radiation 
constant. For the equation of state, the gas pressure is given by the ideal gas law, p = 
RqPT / p, where p is the mean molecular weight and Rq is the gas constant. The temperature 
T is proportional to the specific internal energy, e, by the relation p = (7 — l)pe = RqpT / p, 
where 7 is the specific h eat ratio. To close the system o f equations, we use the fiux-limited 



diffusion approximation (jLevermore fc Pomraning 
-grad Eo, 



19811 ) for the radiative fiux: 



p{k + a 



Fo-- 
and 

Po = Eq ■ TEdd 



(9) 



(10) 



where A and Tgdd are the flux-limiter and the Eddinqt on Tensor, respectively, for which 
we use the approximate formulas given in iKlevi (|l989r ). The formulas fulfill the correct 
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limiting conditions in the optically thick diffusion limit and the optically thin streaming 
limit, respectively. 



For the kinematic viscosity, z/, we adopt a modified version (jPapaloizou fc Stanleylll98fi[ ) 
of the standard a- model. The modified prescription for u is given by 

u = a Cs min [Hp, H] , (11) 

where a is a dimensionless parameter, Cg the local sound speed, and H the disc height. 
Hp = p/ I grad p \ is the mid-plane pressure scale height inside of the disc but the local 
pressure scale height outside of the disc. 



3 NUMERICAL METHODS 



The set of partial differential equations (l)-(6) is numerically solved by a finite-difference 
method under adequate initial and bo undary cond ition s. The numerical sch emes used are 

Kiev (1^3) and 



basically the same as that described by 



Oku da et a1 



( 19971 ). The methods 



are based on an explicit-implicit finite difference scheme. Grid points in the radial direction 
are spaced logarithmically as Ar/r = 0.1, while grid points in the angular direction are 
equally spaced, but more refined near the equatorial plane, typically A( = 7r/150 for 7r/2 ^ 
C ^ 7r/6 and A( = 7r/300 for 7r/6 ^ C ^ 0. Although the radial mesh-sizes do not have a 
fine resolution to examine detailed disc structure, we consider the mesh-sizes to be sufficient 
for examination of the global behaviors of the disc, the jet, and the mass-outflow rate, as is 
later shown in subsection 4.2. 



3.1 Model Parameters 

We consider a Schwarzschild black hole with mass M^, = IOMq and take the inner-boundary 
radius -Rin of the computational domain as 2rg. The model parameters used are listed in table 
1, where -Rmax is the outer boundary radius, and m is the input accretion rate normalized 
to the Eddington critical accretion rate Me which is given by 

Me = 16Le/c2, (12) 

where Le is the Eddington luminosity. L, Mout, and Min in table 1 denote the total luminosity 
through the outer boundary surface, the mass-outflow rate ejected from the outer boundary, 
and the mass-inflow rate swallowed into the black hole through the inner boundary at the 
final state. The viscosity parameter a = 10~^ is adopted, since in the previous examination 
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Table 1. Model parameters 



Model 


rh 


Af(gs-i) 




L/Le 


Mout/M 




1 


3 


8 xlO^^ 


5 X lO* 


variable 


variable 


0.4 


2 


15 


4 xlO^o 


5 X 10** 


8 


0.08 


0.2 



we found that the low value of a = 10 ^ is more advantag eous to 



high velocity jets than the high value case with a = 0.1 (|Okuda 



'ormation of the collimated 



2002). 



3.2 Initial Conditions 

The initial conditions consist of a cold, dense, and optically thick disc and a hot, rarefied, 
and optically thin atmosphere around the disc. The initial disc at r/rg ^ 3 is approximated 
by the Shakura-Sunyaev standard model but, at 2 ^ r/rg ^ 3, it is taken to be sub-Keplerian 
disc with a free-fall velocity at the inner boundary. To examine how the high velocity jet 
evolves through the surrounding matter, it is important, for us, to set up the appropriate 
den sity distribution aro und the disk. Observations and hydrodynamical modeling of SS 433 



jet (jKotani et al 



19961) suggest that the X-ray emitting region of SS 433 would be in the 
range of r = 10^*^ - 10^^ cm from the central source and have gas densities of ~ 10^^^ - 
10^^'^ g cm^^ at the base of the X-ray jets. Focussing on these facts, we take the density at 
the outer boundary radius of 1.5 x 10^^ cm to be ~ 10~^^ g cm~^. Furthermore, assuming 
the density distribution of p ~ r^^, we construct an initial hot rarefied atmosphere around 
the disc to be approximately in radially hydrostatic equilibrium. 



3.3 Boundary Conditions 

Physical variables at the inner boundary, except for the velocities, are given by extrapolation 
of the variables near the boundary. However, we impose limited conditions that the radial 
velocities are given by a free-fall velocity and the angular velocities are zero. On the rotational 
axis and the equatorial plane, the meridional tangential velocity w is zero and all scalar 
variables must be symmetric relative to these axes. The outer boundary at r = -Rmax is 
divided into two parts. One is the disc boundary through which matter is entering from 
the outer disc. At the outer-disc boundary we assume a continuous inflow of matter with a 
constant accretion rate M. The other is the outer boundary region above the accretion disc. 
We impose free-floating conditions on this outer boundary and allow for outflow of matter. 
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Time (Rg/c) 



Figure 1. Time evolution of total luminosity L emitted from the outer boundary surface, disc luminosity from the disc 
surface, mass-outflow rate Mout from the outer boundary, and mass-inflow rate Af;,; swallowed into the black hole through the 
inner boundary for model 1, where time is shown in units of rg/c. 

whereas any inflow is prohibited here. We also assume the outer boundary region above the 
disc to be in the optically-thin limit, |-Fo| ^ cEq. 



4 NUMERICAL RESULTS 
4.1 CaseofM = 3ME 



Model 



model ( Okuda 



with M = 3Me has the same input accretion rate as BH-1 in the previous disc 



20021 1. but the computational domain used is by two orders of magnitude 
larger than the previous one. This is due to the following reasons. We use the initial discs 
based on the Shakura-Sunyaev model. In the super-Eddington accretion regime, the initial 
conditions should be carefully used on treatment of the outer disc-boundary, because we 
don't know a priori reasonable conditions of the very luminous and thick disc. However, if 
the outer boundary radius is taken to be considerably large, this problem will be avoidable 
because the outermost disc is well described by the Shakura-Sunyaev model. Another reason 
is due to a requirement that we want to examine the behavior of high-velocity jet over a 
distant region, which are observed in the X-ray source SS 433. 

The initial disc used is geometrically thick as H/ r ~ 1 in the inner region, but the relative 
disc height H/r and the ratio (3 of the gas pressure to the total pressure are ~ 0.05 - 0.1 
and ~ 0.95 - 0.98, respectively, at the outer disc boundary for models 1 and 2. 
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Figure 2. Velocity vectors and temperature contours in logaritlimic scale on the meridional plane at t = 6 x 10'' (a), 1.5 x 10® 
(b), 2.5 X 10® (c), 2.8 X lO^rg/c (d) for model 1. The contour lines with labels of log T = 5, 6, 7, 8, 9, 10, and 10.5 are shown. 
The reference vector of light is shown by a long arrow in (a). The initially anisotropic high- velocity jet (a) generated in the 
inner region of the disc evolves into the isotropic outflow (d) in the distant region, 'disc' in the figure denotes the cold and 
dense accretion disc. 



Fig. 1 shows the time evolutions of the total luminosity L emitted from the outer bound- 
ary surface, the disc luminosity emitted from the disc surface, the mass-outflow rate Mont 
from the outer boundary, and the mass-inflow rate Min swallowed into the black hole through 
the inner boundary for model 1. Prom the luminosity curve, we can see three characteristic 
stages: stage 1 ( t = - 2.5 x lO^Vg/c ), stage 2 ( t = 2.5 - 5.5 x lO^rg/c ), and stage 3 ( 
t ^ 5.5 X lO^Tg/c ). At the stage 1, matter in the inner region of the disc is ejected strongly 
outward and propagates vertically to the disc plane. The high-velocity jet with a velocity of 
~ 0.2c reaches the outer boundary in the polar direction at i ~ i?niax/0.2c ~ 2.5 x lO^rg/c, 
and then the mass outflow through the outer boundary begins, and the total luminosity L 
attains to ~ 10^° erg s""^. At the stage 2, the luminosity flares up to ~ 4 x 10^° erg s"^ and 
maintains its high value for a time of 3 x lO^rg/c and switches to the stage 3, flaring up by 
more than one order of magnitude. During the stage 3, the luminosity shows quasi-periodic 
variations with a period of 2.5 x lO^rg/c and an amplitude of factor ten. Through these 
stages, the mass-outflow rate Mout becomes as high as ~ 10^^ - 10^^ g s~^, but it is variable 
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Figure 3. Same as figure 2 but at t = 5.1 x 10^ (a), 5.9 x 10^ (b), 8.0 x 10^ (c), 1.1 x 10^ (d) rg/c for model 1. (a) and (c) show 
active outflows in the high temperature state, but (b) and (d) show inflows (no mass-outflow) in the low temperature state. 



and intermittent as is found in Fig.l; that is, the outflow switches occasionally to inflow. 
The mass-inflow rate Min fluctuates rapidly, but the averaged rate is constantly ~ 3 x 10^^ 
g s""^ except the final phases, where becomes very small as ~ 10^^ g s~^. In model 1, the 
heating wave (see Fig. 7) and the density wave, which generate initially in the innermost 
disc region, propagate outward with times and reach the outer disc region of r ~ 3 x lO^rg at 
t ~ 1.1 X lO^rg/c. At these phases, a part of the waves begins to back inward as an ingoing 
wave, and gradually the ingoing density wave leads to lower density in the inner disc region 
and resultantly to smaller mass-inflow rate at the inner boundary. 

L and are given by calculating / FociS, where the surface integral is taken over the 
outer boundary surface and the disc surface, respectively. Although it is difficult for us to 
specify correctly the disc surface in the geometrically thick disc, we define it here as a location 
where the density drops to a tenth of the central density on the equatorial plane. This may 
lead to some errors in estimation of the disc surface and accordingly L^- The complicated 
features of during the stage 3 in Fig.l may be due to unreliable determinations of the 
disc surface. If the atmosphere above the disc is optically thin, must agree with L. 
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has an initial value of ~ 8 x 10^^ erg derived from the Shakura-Sunyaev model, while 
L is initially small and becomes comparable to at the phase t ~ -Rmax/c = 5 x lO'^Vg/c 
when radiation from the disc attains to the outer boundary radius. After t ~ 2.5 x lO^rg/c, 
the disc luminosity arises to ~ 1.3 x 10^° erg s~^ and fluctuates around the value by a 
factor 2 during a time of 2.5 x lO^rg/c, and then increases sharply to 3 x 10^^ erg s~^. The 
total luminosity L also denotes same behavior as the disc luminosity but with a phase lag of 
5 x lO^Tg/c, which is the radiation transit time from the disk surface to the outer boundary 
radius. These time-dependent behaviors of the luminosities are discussed later in terms of 
the thermal instability of the slim disc. 

Fig. 2 denotes the temperature contours with velocity vectors at t = 6 x 10^, 1.5x 10^, 
2.5 X 10^, and 2.8 x lO^rg/c in the stage 1 for model 1. The evolutionary features of outward 
moving jets are found here, where the outermost crowded contours of temperature show the 
shocked region due to the interaction of the jet with the surrounding gas, and the reference 
velocity of hght is indicated by a long arrow in Fig. 2 (a). The temperatures are as high as 
~ 10^ - lO^'^ K in the high-velocity region and 'disc' in the figure denotes the cold and dense 
disc. The initially anisotropic jet (a) propagates vertically to the disc plane but gradually 
loses its anisotropic nature with increasing times, and finally the jet becomes isotropic fiow 
(d) after the jet reaches the outer boundary. The figure also shows that there exist an inactive 
phase (b) and an active phase (c) of outflows. Fig. 3 shows the temperature contours with 
velocity vectors at t ^ 5.0 x 10^ 5.9 x 10^ 8.0 x 10^ and 1.1 x lOVg/c in the stage 3. 
Fig. 3 (a) and (c) show powerful outflows in the high-temperature state above the disc, 
whereas Fig. 3 (b) and (d) denote inflow states in the low-temperature state. The inactive 
inflow states, accompanying with no mass-loss, are likely to occur after a large flare up of 
the luminosity. 

As is found in Figs. 2 and 3, the gas temperatures in the outflow region outside of the 
disc are partly too high as ~ 10^ - 10^^ K. These temperatures may be unrehable because 

we did not take account of other physical processes, such as Compton processes and pair 
production- annihilation, which would be important at such high temperatures. As far as the 
disc is concerned, however, the temperatures are in the range of 10^ ^ T < 10® K. 
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Figure 4. Time evolution of total luminosity L, disc luminosity Lj, mass-outflow rate Mout from the outer boundary, and 
mass-inflow rate Min swallowed into the black hole through the inner boundary for model 2. 





r/Rg 



Figure 5. Velocity vectors and temperature contours in logarithmic scale on the meridional plane at t = 5 x 10* (a), 8 x 10* 
(b), 1.5 X 10"'' (c), 2.5 X 10 ' (d) but for model 2. Differently from model 1, the initially anisotropic high-velocity jet along the 
rotational axis maintains its anisotropic shape when it expands outward. 
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Figure 6. Velocity vectors and contours of density (a) and temperature (b) in tlie inner region of tiie disc and the higli- velocity 
jet at t = 2.5 x 10^ for model 2. The contour labels denote the logarithmic values of the density p (g cm~^) and the temperature 
T(K). 



4.2 Case of M = 15Me 

The time evolutions of L, La, Mout, and Min for model 2 are shown in Fig. 4. At the initial 
phases, the disc luminosity and the mass-inflow rate Min are very dependent on the initial 
structure of the inner disc and they fluctuate largely, but they settle to their steady state 
values in a time-scale of ~ 2 x lO'^Vg/c. The total luminosity L becomes comparable to the 
disc luminosity at t ~ -Rmax/c = 5 x lO'^rg/c when radiation from the disc arrives at the 
outer boundary, and finally it accords with L^. On the other hand, after t ~ i?niax/0.4c = 
lO^rg/ c which is the jet transit time from the disc to the outer boundary, the mass outflow 
begins and Mout gradually settles down to a steady state value of 3 x 10^^ g s~^. Finally, 
L = Ld = 1.3 X 10^° erg s~^. The mass- inflow rate swallowed into the black hole fluctuates 
always by a small factor around 8 x 10^^ g s^^. 

Fig. 5 shows the evolutionary features of the high- velocity jet at t = 5 x 10*^ (a), 8 x 10"^ 
(b), 1.5 X 10^ (c), and 2.5 x 10^ (d). The high-velocity jet propagates vertically to the disc 
plane in the same way as model 1. The jet expands gradually far from the rotational axis 
with increasing times. After the phase in Fig. 5 (b), the jet arrives at the outer boundary 
in the polar direction but the anisotropic nature of the jet is remained even after it passed 
through the outer boundary. If shapes of the jets found in Figs. 5 (a) and (b) are always 
kept throughout the time evolution, we can expect a small coUimation angle of ^ 10° for 
the jet. 
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Fig. 6 shows the velocity vectors and the contours of density p (g cm~'^ ) and temperature 
T (K) of the disc and the high- velocity jet in the inner region at t = 2.5 x lO^rg/c for model 2. 
This shows a rarefied, hot, and optically thin high velocity jet region and a dense, cold, and 
geometrically thick disc region. The high-velocity region and the disc region are bounded 
by the funnel wall, a barrier where the effective potential due to the gravitational potential 
and the centrifugal one vanishes. The jet velocity in the funnel region between the funnel 
wall and the Z-axis amounts to ~ 0.3c. In the inner disc, the flow is convectively unstable. 
In the innermost region of the disc, there exist the inner advection-dominated zone and 
the outer convection-dominated zone. The transition radius rtr between the two zones is 
approximately 10 rg. At r ^ rtr, roughly half the mass at any time at any radius is flowing 
in and flowing out, respectively, and the mass-inflow rate balances the mass-outflow rate. 
Near to the transition region, by way of the convective zones and the accretion zone, matter 
is accreted towards the equatorial plane. The accreting matter, which is carried to the 
transition by convection, partly diverts into the high-velocity region and partly flows into 
the inner advection-dominated zone. The gas diverged into the optically thin high-velocity 
region is originally subsonic but is soon accelerated by the strong radiation field up to a 
relativistic velocity. These features of the accreting and outflow materials in the inner region 
are also found in model 1, even when it is at an inflow phase in the outer region. 

In Fig. 6, we see that at the disc boundary the densities lie between label — 6 and — 
7 and a temperature of ~ 10^ K in the disc jumps to ~ 10^ K. The opening angle of the 
high- velocity jet to the Z-axis is ~ 30°. These density and temperature profiles, the jet 
velocity, and the ope ning angle sh ow good agreements with those in the previous calculation 



extending to 400r„ (jOkuda 



20021 ) ■ where the relative radial mesh of Ar/r = 0.03 is finer 



than 0.1 in the present one but the mesh-sizes in the angular direction have same resolution. 
Therefore, the present results over an extended region of ~ 5 x lO^'rg are considered to be 
also valid in spite of the coarse radial mesh-sizes. 

The disc and the jet in model 2 attain almost to their steady states, while model 1 shows 
variable luminosities and non-steady gas flows. These behaviors of the disc luminosity in 
models 1 and 2 are considered in terms of the thermal instability of the slim accretion disc 
model. The thermal instability of the standard disc is interpreted by the relation between the 
accretion rate m and the surface density S at a fixed radius r. When the curve of m = m(E) 
has a characteristic S-shape with three branches (lower, middle, and upper), the discs in the 
upper and lower branches are stable against thermal instability, but in the middle branch 
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Rg 



Figure 7. Temperature versus radius on the equatorial plane at t= (line), 2.5 X 10^ (dashed line), 5.1 X 10^ (dotted line), 
8.0 X 10^ (dash-dot line), and 1.1 X lO" (dash-dot-dot line) Tg/c for model 1. 




r/Rg 



Figure 8. Disc height H versus radius on the equatorial plane at t= (line), 2.5 X 10^ (dashed line), 5.1 X 10^ (dotted line), 
8.0 X 10^ (dash-dot line), and 1.1 X 10^ (dash-dot-dot line) rg/c for model 1. 



it is unstable. If the disc is under the unstable middle branch, it is expected that the disc 
exhibits a limit-cycle behavior of lumin osity. 



Taking account of advective cooling, 



Abramowicz et al, 



( 19881 ) found the S-shaped m — S 



curve at r = 5rg, corresponding to moderately super-Eddington accretion rates from the 
transonic solutions of a black hole with = IOMq and a = 10^^. According to their results, 
the inner disc with m = 15 in our model 2 corresponds to the stable upper branch, while 
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the disc with m = 3 belongs marginally to the unstable middle branch and it locates near 
the turning point between the upper and middle branches. Fig. 7 denotes the temperature 
versus radius on the equatorial plane at t = 0, 2.5 x 10^, 5.1 x 10^, 8.0 x 10^, and 1.1 x 



lO^rg/c for model 1. The maximum luminosity appears at t = 5.1 x 10^rg./c. The evo 



of the disc temperatu r e is q ualitatively similar to the results by 



Honma et al 



ution 



(19911) and 



Szuszkiewicz fc Miller! fll99a l. who found a limit-cycle behavior of the slim disc with m 



0.06 and a = 0.1 by one-dimensional calculations of vertically integrated equations. The 
model with m = 0.06 and a = 0.1, which belongs to the unstable middle branch of m - 
S curve at r/rg = 5, locates just near another turning point between the unstable middle 
branch and stable lower branches. In their results, the outgoing heating wave reaches about 
lOOrg, while in ours the heating wave attains to ~ 3 x lO^rg. The disc height H versus 
radius at same phases as Fig. 7 is shown in Fig. 8. In the initial phase, H/r is ~ 0.1 at 
r/rg ^ 10'^, but the outgoing heating wave leads to H/r ~ 1 in the o uter disc region and 



the radiation pressure is do r ninan t there. Compared with the results by 



and 



Honma et al, 



(|l99li l 

Szuszkiewicz fc Milled ()l998( ). the heating wave reaches the extensively distant region 



because of the much higher input accretion rate used. Apart from the effects due to the 
much higher accretion rate, we can see similar characteristics of the disc evolution between 
theirs and us and expect a limit cycle behavior of the luminosity in model 1. However, in 
order to confirm the limit cycle phenomenon of the luminosity, the computational time in 
model 1 may not be sufficient because the temperature evolution shows only the phases of 
the outgoing heating wave. 

On the other hand, the super-Eddington disc with m = 15 is stable against the thermal 
instability due to the strong advective cooling. In this model, the outgoing heating wave with 
temperatures of 10^ K ^ T ^ 10^° K generates near the inner edge of the disc at the initial 
phases and reach about ~ lO^rg. However, these waves fade out soon and successively the 
second heating wave generates near the inner edge and oscillates at 2 ^ r /rg ~ 10 during a 
time of ~ lO^rg/c, and finally they are stabilized completely. 



5 COMPARISON WITH SS 433 

SS 433 is a typical stellar black-hole candidate at a highly super-Eddington luminosity. Our 
results may be compared with the observations in SS 433. The observations show many 
emission lines of heavy elements, such as Fe and Ni, which denote the unique red and blue 
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Doppler shifts of 0.26 c in the X-ray region. The colhmation angle of SS 433 jets seems to be 
as sr nall as 0.1 radian (sev eral degrees). From an observational constraint of the mass-outflow 
( 19961 ) shows Mout ^ 3 x 10^^ g s'^ in SS 433. 



rate. 



Kotani et al 



The unique velocity, 0.26 c, may be reasonably explained in terms of the relativistic ve- 
locities accelerated by the radiation-pressure force in the inner region of the super-Eddington 
disc. In our models, the collimation degree of the high velocity jets will be dependent on the 
initial density distribution of the matter around the disc. However, the density distribution 
used is not a special type of profile but a usual spherical one with p oc r~^. In model 1, 
the initial jet is strongly collimated in the polar region, but it becomes isotropic flow in the 
distant region of ~ 10^^ cm, where the X-ray jet is observed in SS 433. While in model 2 
the small collimation angle of the jet is maintained throughout the jet propagation. If the 
anisotropic nature of the jet found in Fig. 5 (a) and (b) is maintained throughout the disc 
evolution, it is expected that the collimation angle of the jet is as small as ^ 10°. Fur- 
thermore, in addition to models 1 and 2, we examined the intermediate case with m = 10 
between them and found that the case also shows anisotropic outflows same as model 2 
but with a smaller collimation angle than that in model 2. From these results, we speculate 
that the collimation angle may be smaller if we consider other models with much higher 
input-accretion rates. This should be confirmed in future. Of course, axisymmetric MHD 
disc simulations have been used successfully i n the past to reproduce jets formation and 



collimation of the jets since 



Shibata &: Uchidal (|l986[ ). and the recent 3-D MHD simulations. 



which have directed at the internal dynamics of the disc itself and the resulting accretion 
flow, also show an unbound outflowing jet along the Z axis confined by magnetic pressure. 
However, as far as the very luminous X-ray source SS 433 is concerned, the jet acceleration 
is considered to be mainly in thermal origin but not nonthermal one. Therefore it is very 
advantageous for us to be able to reproduce a small collimation angle of the jet in terms of 
the thermal model considered here. 



The mass-outflow rate in model 1 is as large as ~ 10^^ - 10^^ g s^^ but is variable and 
intermittent with times, while in model 2 it is as large as ~ 3 x 10^^ g s~^, and it is stable 
and comparable to the mass-outflow rate estimated in SS 433. 
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6 CONCLUDING REMARKS 



Black-hole accretion discs and jets with 3 and 15 Me have been examined by time- dependent 
two-dimensional radiation hydrodynamical calculations over a far distant region and a long 
time-scale. Starting with the initial disc based on the Shakura-Sunyaev disc model, we ex- 
amined the time evolutions of the disc and the jet. As the result, the inner disc becomes 
geometrically thick, optically thin, and convectively unstable. The dominant radiation pres- 
sure force in the inner region of the disc accelerates the gas flow vertically to the disc plane, 
and the high-velocity jets with 0.2 - 0.4c are formed along the rotational axis. In model 1 
with the lower accretion rate, the initially anisotropic high-velocity jet becomes gradually 
isotropic flow in the distant region of ~ 10^^ cm. The mass outflow occurs through the 
entire outer boundary and the mass-outflow rate is large as ~ 10^^ - 10^^ g s~^, but it is 
variable and intermittent with times; that is, the outflow switches occasionally to inflow. At 
the initial phase the luminosity fluctuates around ~ lO^'' erg s~^ by a factor 2 and then it 
flares up to 4 x 10^° erg s~^ and again to ~ 10^^ erg s~^ after a time of 2.5 x lO^rg/c. After 
then the luminosity behaves like a periodic variation with a period of ~ 2.5 x lO^rg/c and 
an amplitude of factor ten. 

On the other hand, the jet in model 2 with the higher accretion rate maintains its initial 
anisotropic shape even after the jets go far away, and the mass outflow occurs only through 
the polar direction. The disc luminosity and the mass-outflow rate attain to their steady 
values of 1.3 x 10^° erg s~^ and 3 x 10^^ g s~^, respectively. These different behaviors of the 
unstable disc luminosity in model 1 and the stable disc luminosity in model 2 are interpreted 
in terms of the thermal instability of the slim accretion disc model. The super- Eddington 
model with 15 Me is promising to explain a small collimation degree of the jet and a large 
mass-outflow rate observed in the X-ray source SS 433. 

It shoud be noticed that these conclusions are dependent on the viscosity parameter 
a because we examined only the low-viscosity flows with a = 10~^. In the previous nu- 
merical simulations, we found that the low-viscosity flow is more advantageous to form a 
well coUimated powerful jet than the high- viscosity flow with a = 0.1. The characteristics 
of the accretion flows for different a can qualitatively explained by the role of viscosity in 



of low viscosity ( a ^ 0.01) are generally convection-dominated flows ( 


Stone 


et al. 


1999; 


leumenshchev & Abramowicz 


1999; 


Naravan et al. 


2000; 


Ouataert & Gruzinov 


2000 


). The 
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convective instability results in the vortex motion of gas of different spatial scales. Con- 
vection cells stretch from the disc mid-plane to the disc surface in the inner region of the 
disc. The convective flow is subsonic in the mid-plane but becomes comparable to the sonic 
speed near to the disc surface. As a result, a part of the convective flow escapes from the 
disc as the disc wind, dissipating energy of powerful accoustic waves produced by convection 
in the disc. The outflowing material originates not only from the innermost disc near the 
transition region between the inner advection-dominated and outer convection-dominated 
zones but also from the outer upper disc as the disc wind. The disc wind is incorporated 
into the rarefied, hot, and high- velocity region and is acclerated due to the strong radiation 
field in the rarefied region. However, in the case of a high viscosity, the strong accretion flow 
smoothes the convective instability and is globally stable to the instabihty. As the result, 
the powerful jet as is found in the low- viscosity case is not reproduced. At present, we have 
no clear conclusion as to which viscosity case to be realized in actual black hole accretion 
flows. Future 3-D MHD simulations of the black hole accretion flows may give a decisive 
answer to the viscosity problem. 
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